Quantum information processing: 
The case of vanishing interaction energy 

Miroljub Dugic 1 '* and Milan M. Cirkovic 2 '** 

* Department of Physics, Faculty of Science 
P. 0. Box 60, 34000 Kragujevac, Serbia 
** Astronomical Observatory, Volgina 7, Belgrade, Serbia 
1 E-mail: dugic@knez.uis. kg. ac.yu 
2 E-mail: mcirkovic@aob.aob. be;. ac.yu 

O ■ 

Abstract: We investigate the rate of operation of quantum "black boxes" ("oracles") 
and point out the possibility of performing an operation by a quantum "oracle" whose 
average energy equals zero. This counterintuitive result not only presents a generalization 
of the recent results of Margolus and Levitin, but might also sharpen the conceptual 
r**- ! distinction between the "classical" and the "quantum" information. 

^ ! PACS: 03.67.Lx, 03.65. Ud, 03.65.Yz 

1. Introduction: information-processing bounds 

o : 

In the realm of computation, one of the central questions is "what limits the laws 
O ■ of physics place on the power of computers?" [1,2]. The question is of great relevance to 
a wide range of scientific disciplines, from cosmology and nascent discipline of physical 
Qh! eschatology [3-5] to biophysics and cognitive sciences which study information processing 
in the conscious mind [6, 7]. One of the physical aspects of this question refers to the 
minimum time needed for execution of the logical operations, i.e. to the maximum rate 
of transformation of state of a physical system implementing the operation. From the 
fundamental point of view, this question tackles the yet-to-be-understood relation between 
• £h ! the energy (of a system implementing the computation) on the one, and the concept of 
information, on the other side. Apart from rather obvious practical interest stemming 
from the explosive development of information technologies (expressed, for instance, in 
the celebrated Moore's law), this trend of merging physics and information theory seems 
bound to offer new insights into the traditional puzzles of physics. Specifically, answering 
the question above might shed new light on the standard "paradoxes" of the quantum 
world [8-10]. 

Of special interest are the rates of the reversible operations (i.e. of the reversible 
quantum state transformations). To this end, the two bounds for the so-called "orthogonal 
transformations (OT)" are known; by OT we mean a transformation of a (initial) state 

to a (final) state while (^i\^f) = 0. First, the minimum time needed for OT 

can be characterized in terms of the spread in energy, AH, of the system implementing the 
transformation [11-15]. However, recently, Margolus and Levitin [16, 17] have extended 
this result to show that a quantum system with average energy (H) takes time at least 
t = h/4((H) — Eq) to evolve to an orthogonal state, where Eq is the ground state energy. 
In a sense, the second bound is more restrictive: a system with the zero energy (i.e. in the 
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ground state) cannot perform a computation ever. This however stems nothing about the 
nonorthogonal evolution which is still of interest in quantum computation. 

Actually, most of the efficient quantum algorithms [f 8-21] employ the so-called quan- 
tum "oracles" (quantum "black boxes") not requiring orthogonality of the initial and the 
final states of the composite quantum system "input register + output register (I + O)" 
[22-24]. Rather, orthogonality of the final states of the subsystems' (e.g. O's) states is 
required, thus emphasizing a need for a new bound for the operation considered. 

In this paper we show that the relative maximum of the rate of operation of the 
quantum "oracles" may point out the zero average energy of interaction in the composite 
system I + O. Actually, it appears that the rate of an operation cannot be characterized 
in terms of the average energy of the composite system as a whole. Rather, it can be 
characterized in terms of the average energy of interaction Hamiltonian. Interestingly 
enough, the ground state energy Eq plays no role, here, and the absolute value of the 
average energy of interaction (\(Hi nt )\) plays the role analogous to the role of the difference 
(Hint) —Eq in the considerations of OT. Physically, we obtain: the lower the average energy, 
the higher the rate of operation. This result is in obvious contradistinction with the result 
of Margolus and Levitin [16, 17] — in terms of the Margolus-Levitin theorem, our result 
would read: the lower the difference (H- int ) — Eq, the higher the rate of (nonorthogonal) 
transformation. On the other side, our result is not reducible to the previously obtained 
bound characterized in terms of the spread in energy [11-15], thus providing us with a new 
bound in the quantum information theory. 

2. The quantum "oracle" operation 

We concern ourselves with the bounds characterizing the rate of (or, equivalently, the 
minimum time needed for) the reversible transformations of a quantum system's states. 
Therefore, the bounds known for the irreversible transformations are of no use here. Still, 
it is a plausible statement that the information processing should be faster for a system 
with higher (average) energy, even if — as it is the case in the classical reversible information 
processing — the system does not dissipate energy (e.g. [25]). This intuition of the classical 
information theory is justified by the Margolus-Levitin bound [16, 17]. However, this 
bound refers to OT, and does not necessarily applies to the nonorthogonal evolution. 

The typical nonorthogonal transformations in the quantum computing theory are 
operations of the quantum "oracles" employing "quantum entanglement" [18-21, 23, 24]. 
Actually, the operation considered is defined by the following state transformation: 

\*i)io = C *\ x )i ® -> \*f)io = ]TC*|£>7 ® \f(x)) , (1) 

X X 

where {!#)/} represents the "computational basis" of the input register, while \0)o repre- 
sents an initial state of the output register; by " /" we denote the oracle transformation. 

The key point is that the transformation (1) does not [18-21] require the orthogonality 
condition io(^i\^f) io = to be fulfilled. Rather, orthogonality for the subsystem's states 
is required [18-21]: 

o{f(x)\f(x')) o = 0, x^x' (2) 
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for at least some pairs (x, x'), which, in turn, is neither necessary nor a sufficient condition 
for the orthogonality io{^i\^ f)io = to be fulfilled. 

Physical implementation of the quantum oracles of the type Eq. (1) is an open ques- 
tion of the quantum computation theory. However (and in analogy with the quantum 
measurement and the decoherence process [26-31]), it is well understood that the imple- 
mentation should rely on (at least indirect, or externally controlled) interaction in the 
system / + O as presented by the following equality: 

\* f ) IO = U(t)\*i) IO = U(t)J2C x \x)i ® \0) o = C *\ x )i ® *))o> ( 3 ) 

X X 

where U (t) represents the unitary operator of evolution in time (Schrodinger equation) for 
the combined system 7 + 0; index t represents an instant of time. Therefore, the operation 
(1) requires the orthogonality: 

o{f(x,t)\f(x',t)) o = 0, (4) 

which substitutes the equality (2). 

So, our task reads: by the use of Eq. (4) , we investigate the minimum time needed for 
establishing of the entanglement present on the r.h.s. of Eq. (1), i.e. of Eq. (3). 

3. The optimal bound for the quantum oracle operation 

In this Section we derive the optimal bound for the minimum time needed for execution 
of the transformation (1), i.e. (3), as distinguished by the expression (4). We consider 
the composite system "input register + output register (I + 0) v defined by the effective 
Hamiltonian: 

H = H I +H Q +H int (5) 

where the last term on the r.h.s. of (5) represents the interaction Hamiltonian. For 
simplicity, we introduce the following assumptions: (i) dH/dt = 0, (ii) [Hi,Hi nt ] = 0, 
[Ho, Hi nt ] = 0, and (iii) H int = CAj <S>Bo, where Aj and Bo represent unspecified observ- 
ables of the input and of the output register, respectively, while the constant C represents 
the coupling constant. As elaborated in Appendix I, the simplifications (i)-(iii) are not very 
restrictive. For instance, concerning point (i) — widely used in the decoherence theory — 
one can naturally relax this condition to account for the wide class of the time dependent 
Hamiltonians; cf. Eq. (1.1) of Appendix I. Another way of relaxing the condition (i) is to 
assume the sudden switching the interaction in the system on and off. 

3.1 The bound derivation 

Given the above simplifications (i)-(iii), the unitary operator U (t) (cf. Eq. (3)) spectral 
form reads: 

U(t) = ^ exp{-zt(e x + Ei + C lxi )/fi}P Ix ® ft 0i . (6) 

x,i 
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The quantities in Eq. (6) are defined by the following spectral forms: Hi = J2 X e x Pi Xl 

Ho = J2i EiKoi, and H 1Qt = CJ2 X:i lxiPi x <S> n 0i ; bearing in mind that Aj = J2 X a x Pi x 

and B = Y.i b i^Oi, the eigenvalues j xi = a x bi. 

^From now on, we take the system's zero of energy at the ground state by the exchange 
E x i — > E x i — E Q ; E x i = e x + Ei + C^ x i, E Q is the minimum energy of the composite 
system — which Margolus and Levitin [16, 17], as well as Lloyd [1], have used. Then one 
easily obtains for the output- register's states: 

|/0r, t))o = J2 ex P{-^( e - + E i + Clxi - E o )/ft}n o *|0> o . (7) 

i 

Substitution of Eq. (7) into Eq. (4) directly gives: 

D xx >(t) = (f(x,t)\f(x',t)) = exp{-zt(e x -e x ,)/h}x 
x ^2piexp{-iCt(a x - a x ')bi/h} = 0, = 1, (8) 

i i 

where pi =o (0\Hoi\0)o- The expression (8) represents the condition for "orthogonal 
evolution" of subsystem's (O's) states bearing explicit time dependence; the ground energy 
E does not appear in (8). 

But this expression is already known from, e.g., the decoherence theory [26-31]. Ac- 
tually, one may write: 

D X x'{t) = exp{-it(e x - e x >)/h}z xx '(t), (9) 

where 

z X x'(t) = ^2piexp{-iCt(a x - a x >)bi/h} (10) 

i 

represents the so-called "correlation amplitude", which appears in the off-diagonal elements 
of the (sub)system's (O's) density matrix [26]: 

POxx'(t) = C x C*,Z xx '(t). 

So, we could make a direct application of the general results of the decoherence theory. 
However, our aim is to estimate the minimum time for which D xx i (t) may approach zero, 
rather than calling for the qualitative limit of the decoherence theory [26]: 

lim \z xx ,(t)\ = 0, (11) 

t— >oo 

or equivalently lim^oo z xx i (t) — > 0. 

In order to obtain the more elaborate quantitative results, we shall use the inequality 
cos a; > 1 — (2/n)(x + sinx), valid only for x > [16, 17]. However, the use cannot be 
straightforward. 
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Namely, the exponent in the "correlation amplitude" is proportional to: 

{a x -a x ,)b h (12) 

which need not be strictly positive. That is, for a fixed term a x — a x > > 0, the expression 
Eq. (12) can be both positive or negative, depending on the eigenvalues b^. For this reason, 
we will refer to the general case of eigenvalues of the observable Bo, {bi, — Pj}, where both 
b h (3j > 0,Vi, j. 

In general, Eq. (10) reads: 

***'(0 = 42(0 + 42(*). (13a) 

where 

z xx' = ^P;exp{-zCt(a x - a x ,)bi/h}, (136) 

i 

z xx> = J^Pj exp{zCt(a x - a x >)f3j/h}, (13c) 
j 

while J2iPi + Y^jP'j = !• Since botn ( a x - a x ')bi > 0, (a x - a x >)f3j > (Vz, j), one may 
apply the inequality mentioned above. 

Relaxed equality (4) — or relaxed equality (11) — is equivalent to Rez xx ' = and 
lmz xx ' = 0. Now, from Eq. (13a-c) it directly follows: 

Rez xx > = cos[C(a x - a x >)bit/h] 



+ ^p' j cos[C{a x -a x ,)p j t/h], (14) 
j 

which, after applying the above inequality gives: 

4 2 

Rez xx / > 1 - -C(a x - a x ,)(Bi + B 2 )t Imz xx /- 

n tv 

--y^PiSm[C(a x -a x >)bit/h], (15) 

TV z — ' 

i 

where B x = ^pA, and B 2 = Y.jP'jPj- 

Since | J2 i p i sm[C(a x — a x >)bit/h]\ < YliiPi = a < 1, V£, from Eqs. (11) and (15) it 
follows: 

2 4a 4 
Re^ xx / + -Tmz xx > > 1 T C(a x - a x >)(5i + B 2 )t. (16) 

TV TV II 

/.From (16) it is obvious that the condition imposed by Eq. (4) cannot be fulfilled in 
time intervals shorter than T xx r. 

(1-Aa/Tv)h 

XX 4C(a x -a x ,)(B 1 + B 2 y (i7j 
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The expression is strictly positive for a < 7r/4, and which directly defines the optimal 
bound Tent as: 

Tent = SUp {r xx /}. (18) 

The assumption a < tt/4 is not very restrictive. Actually, above, we have supposed 
that neither J2iPi — 1) nor YljPj — 1> wnue the former is automatically satisfied with the 
condition a < n/4. 

3.2 Analysis of the results 

The bound r en t is obviously determined by the minimum difference a x — a x >. This 
difference is virtually irrelevant (in the quantum computation models, it is typically of the 
order of K). So, one may note, that the bound in Eq. (18) may be operationally decreased 
by the increase in the coupling constant C and/or by the increase in the sum B\ + B 2 . As 
to the former, for certain quantum "hardware" [32, 33], the coupling constant C may be 
partially manipulated by experimenter. On the other side, similarly — as it directly follows 
from the above definitions of B\ and B 2 — by the choice of the initial state of the output 
register, one could eventually increase the rate of the operation by increasing the sum 
Si + B 2 . 

Bearing in mind the obvious equality: 

(H int ) = C{A I ){B ) = C(A I )(B 1 - B 2 ), (19) 

one directly concludes that adding energy to the composite system as a whole, does not 
necessarily increase the rate of the operation considered. Rather, the rate of the operation 
is determined by the absolute value of the average energy of interaction, |(Hj nt )|. F° r 
instance, if B x ^ while B 2 = (or B 2 ^ 0, Hi = 0), from Eq. (19) it follows that 
the increase in Bi (or in H 2 , and/or in the coupling constant C) coincides (for (A/) ^ 0) 
with the increase in |(if int )|, as well as with the decrease in the bound Eq. (18). This 
observation is in accordance with the Margolus-Levitin bound [16, 17]: the increase in the 
average energy (of interaction) gives rise to the increase in the rate of the operation (still, 
without any restrictions posed by the minimum energy of either the total, or the interaction 
Hamiltonian) . Therefore, the absolute value |(i?mt)| pl^ys, in our considerations, the role 
analogous to the role of the difference (Hint) — E in the considerations of the "orthogonal 
transformations" . 

However, for the general initial state of the output register, both Bi ^ and B 2 ^ 0. 
Then, e.g., for Bi > B 2 : 

Bt + B 2 = Si(l + «) < 2B U k < 1, (20) 

which obviously determines the relative maximum of the rate of the operation by the 
following equality: 

B 1 = B 2 , k = 1, (21a) 
which, in turn (for (Aj) ^ 0), is equivalent with: 

(H in t) = 0. (216) 
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But this result is in obvious contradistinction with the result of Margolus and Levitin 
[16, 17]. Actually, the expressions (21a,b) stem that, apart from the concrete values of 
Bi and B 2 , the relative maximum of the rate of the operation requires (mathematically: 
implies) the zero average energy of interaction, (Hint) = — which (as distinguished above) 
is analogous to the equality (Hint) — E = for "orthogonal transformations." 

4. Discussion 

Intuitively, the speed of change of a system's state should be directly proportional 
to the average energy of the system. This intuition is directly justified for the quantum 
OT by the Margolus-Levitin theorem [16, 17]. Naively, one would expect this statement 
to be of relevance also for the nonorthogonal evolution. Actually, in the course of the 
orthogonal evolution, the system's state "passes" through a "set" of nonorthogonal states, 
thus making nonorthogonal evolution faster than the orthogonal evolution itself. (At least 
this physical picture is justified for "realistic" interpretations of quantum mechanics, like 
the dynamical-reduction or many- worlds theories.) 

This intuition is obviously incorrect for the cases studied. In a sense, the expressions 
(21) state the opposite: the lower the difference Bi — B 2 (i.e. the lower the average energy 
of interaction, |(Hi n t)|), the faster the operation considered. Therefore, our main result, 
Eq. (21), is in obvious contradistinction with the conclusion drawn from the Margolus- 
Levitin bound [16, 17]: the zero energy quantum information processing is possible and, in 
the sense of Eq. (21), is even preferable. From the operational point of view, the bound 
Tent can be decreased by manipulating the interaction in the combined system I + O, as 
well as by the proper local operations (e.g., the proper state preparations increasing the 
sum Bi + B2) performed on the output register. 

As it can be easily shown, the increase in the sum Bi + B 2 coincides with the increase 
in the spread in Bo, A£?o, i-e. with the increase in the spread AHi nt . This observation 
however cannot be interpreted as to suggest reducibility of the bound in Eq. (17) onto 
the bound characterized in terms of the spread in energy [11-15] — in the case studied, 
AH 1Ilt . Actually, as it is rather obvious, the increase in the spread AH[ nt does not pose 
any restrictions on the average value (Hint)- Therefore, albeit having a common element 
with the previously obtained bound [11-15], the bound in Eqs. (17) and (18), represents a 
new bound in the quantum information theory.* 

From Eq. (17), one directly determines the absolute maximum of the rate of the 
operation, i.e. the absolute minimum of the r.h.s. of Eq. (17). Actually, for Bo bounded 
(which is generally the case for quantum computation models), the inequality B m i n < 
B\ + B2 < -Bmax determines the absolute minimum of the r.h.s. of Eq. (17): 

(l-4a/7r)h 
4C(a x — a x ')B max ' 

* This bound is of interest also for the decoherence theory, but it does not provide us 
with magnitude of the "decoherence time" , td- Actually, one may write — in our notation — 
that td oc (a x — a x >)~ 2 , while — cf. Eq. (17) — r ent oc (a x — a^/) -1 , which therefore indicates 
td 3> Tent- This relation is in accordance with the general results of the decoherence 
theory: the entanglement formation should precede the decoherence effect. 
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where B max (B m i Q ) is the maximum (minimum) in the set {bi,fij}. Interestingly enough, 
the minimum Eq. (22) is achievable also in the following special case: if B\ + B<i = 
Pi&max +Pi0max, while pi = pi = 1/2 and 6 max = /3 max = B max , one obtains, again, that 
(Hint) = 0; by 6 max (Anax) we denote the maximum in the set {bi} ({/%})• 

It cannot be overemphasized: the zero (average) energy quantum information pro- 
cessing is in principle possible. Moreover, the condition (Hint) = determines the relative 
maximum of the operation considered. But this result challenges our classical intuition, 
because it is commonly believed that the efficient information processing presumes an 
"energy cost". In the classical domain, this was established in 1960s by Brillouin [25], 
following the ground-breaking studies of Szilard and others on the problem of Maxwell's 
demon. So, one may wonder if " saving energy" might allow the efficient information pro- 
cessing ever. Without ambition to give a definite answer to this question, we want to 
stress: as long as the "energy cost" in the classical information processing (including the 
quantum- mechanical "orthogonal evolution") is surely necessary, this need not be the case 
with the quantum information processing, such as the entanglement establishing. Actually, 
the "classical information" refers to the orthogonal (mutually distinguishable) states, while 
dealing exclusively with the orthogonal states is a basis of the classical information pro- 
cessing [23]. Nonorthogonal states (i.e. nonorthogonal transformations) we are concerned 
with necessarily refer to the nonclassical information processing. So, without further ado, 
we stress that Eq. (21) exhibits a peculiar aspect of the quantum information (here: of the 
entanglement formation), so pointing to the necessity of its closer further study. 

The roles of the two registers (/ and O) are by definition asymmetric, as obvious from 
Eqs. (1) and (3). This asymmetry is apparent also in the bound given in Eq. (17), which 
is the reason we do not discuss in detail the role of the average value (Aj). Having in 
mind the considerations of Section 3, this discussion is really an easy task not significantly 
changing the conclusions above. 

Finally, applicability of the bound (17) for the general purposes of the quantum com- 
puting theory is limited by the defining expression Eq. (1). The bound in Eq. (17) is of 
no use for the algorithms not employing quantum entanglement. As such an example, we 
may consider Grover's algorithm [34], which does not employ quantum entanglement in 
its oracle operation. As another example, we mention the so-called "adiabatic quantum 
computation" (AQC) model [35, 36]. This new computation model does not employ any 
"oracles" whatsoever. Moreover, the AQC algorithms typically involve non-persistent en- 
tanglement (as distinct from those in Eq. (1)) of states of neighbor qubits (cf. Eq. (II. 1) 
in Appendix II). Therefore, the bound in Eq. (17) is of no direct use in AQC, and cannot 
be used for analyzing this non-circuit model of quantum computation (cf. Appendix II). 

The work on application of Eq. (17) in optimizing entangling circuits is in progress, 
and will be published elsewhere. 

5. Conclusion 

We show that the zero average energy quantum information processing is theoretically 
possible. Specifically, we show that the entanglement establishing in the course of oper- 
ation of some typical quantum oracles employed in the quantum computation algorithms 
distinguishes the zero average energy of interaction in the composite system " input register 
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+ output register" . This result challenges our classical intuition, which plausibly stems a 
need for the "energy cost" in the information processing. To this end, our result, which 
sets a new bound for the nonorthogonal evolution in the quantum information processing, 
sets a new quantitative relation between the concept of information on the one, and of the 
physical concept of energy, on the other side — the relation yet to be properly understood. 
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Appendix I 

Relaxing the simplifications (i)-(iii) of Section 3 does not lead to significant changes 
of our results. This can be seen by employing the arguments of Dugic [29, 30], but for 
completeness, we briefly outline the main points in this regard. 

First, for a time dependent Hamiltonian, which is still a "nondemolition observable", 
[Hit), H(t')] = 0, the spectral form [30]: 



H = ^ lxi (t)P Ix ®iloi. (1-1) 



This is a straightforward generalization of the cases studied and also a wide class of the time 

t 

dependent Hamiltonians. E.g., from (1. 1) it easily follows that the term a x i(t) = J , y x i(t')dt' 

o 

substitutes the term ^ xi t in the exponent of the expression (6) . Needless to say, this relaxes 
the constraint (i) of Section 3 and makes the link to the realistic models of the quantum 
" hardware" . 

To this end, it is worth emphasizing: in the realistic models one assumes the actions 
performed on the qubits in order to design the system dynamics. Interestingly enough, 
such actions usually result in the (effective) time independent model Hamiltonians [23, 
33]. Moreover, some time dependent models allow direct applicability of the above notion; 

e.g., for the controlled Heisenberg interaction, H = J(t)S\ ■ S2, the action reads: U(t) = 

t 

exp(—iK(t)Si ■ S2), where K{t) = f J{t')dt' . As an illustration of the models not fitting 



(i), one can consider NMR models which, in turn, are known to be of only limited use in 
the large-scale quantum computation [23, 32, 33]. We conclude that the realistic models of 
the large-scale quantum computation fit with the relaxed point (i) of our considerations. 

Similarly, relaxing the exact compatibilities (cf. the point (ii) in Section 3) leads to the 
approximate separability — i.e., in Eq. (1.1), there appear terms of small norm — which does 
not change the results concerning the " correlation amplitude" z xx i (t) [26] , and consequently 
concerning D xx / (t) . 

Finally, generalization of the form of the interaction Hamiltonian (cf. point (iii) of 
Section 3) does not produce any particular problems, as long as the Hamiltonian is of (at 
least approximately) separable kind, and also a "nondemolition observable". E.g., from 
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the general form for H int [30], H\ nt = J2 k CkAj}~®Bok, one obtains the term J2k Cfc( a fcx — 
akx')bki, instead of the term of Eq. (12). 

The changes of the results may occur [30] if the Hamiltonian of the composite system 
is not of the separable kind and/or not a " nondemolition observable"; for an example see 
Appendix II. 

For completeness, we notice: a composite-system observable is of the separable kind if 
it can be proved diagonalizable in a noncorrelated (the tensor-product) orthonormal basis 
of the Hilbert space of the composite system [30]. 

Appendix II 

By " nonpersistent entanglement" we mean the states of a composite system which 
can be written as: 

m = J2Cit\it)\it), (II.l) 

i 

i.e. states whose Schmidt (canonical) form is labeled by an instant of time, t (continuously 
varying with time). The occurrence of such forms for AQC can be easily proved by the 
use of the method developed in Ref. [30] applied to, e.g., Eq. (3.5) of Ref. [35]. Needless 
to say, states of the Eq.(II.l)-form are exactly what should be avoided in the situations 
described by Eq. (1). 

To this end, the problem addressed in Ref. [30] reads: "what characteristics of the 
system Hamiltonian are required in order to attain the persistent entanglement (cf. Eq. 
(1))?". The answer is given by the points (i)-(iii) (but see Appendix I). In other words, 
as long as the conditions (i)-(iii) are fulfilled, nonpersistent entanglements do not occur 
in the system. As a corollary, having (i)-(iii) in mind, the nonpersistent entanglement of 
AQC cannot be (at least not directly) addressed within the present considerations. 
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